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Myosin V (myoV) is a two-headed myosin capable of taking many successive steps along actin per diffusional en- 
counter, enabling it to transport vesicular and ribonucleoprotein cargos in the dense and complex environment 
within cells. To better understand how myoV navigates along actin, we used polarized total internal reflection fluo- 
rescence microscopy to examine angular changes of bifunctional rhodamine probes on the lever arms of single 
myoV molecules in vitro. With a newly developed analysis technique, the rotational motions of the lever arm and 
the local orientation of each probe relative to the lever arm were estimated from the probe's measured orientation. 
This type of analysis could be applied to similar studies on other motor proteins, as well as other proteins with 
domains that undergo significant rotational motions. The experiments were performed on recombinant con- 
structs of myoV that had either the native-length (six IQ, motifs and calmodulins [CaMs]) or truncated (four IQ 
motifs and CaMs) lever arms. Native-length myoV-6IQ mainly took straight steps along actin, with occasional small 
azimuthal tilts around the actin filament. Truncated myoV-4IQ showed an increased frequency of azimuthal steps, 
but the magnitudes of these steps were nearly identical to those of myoV-6IQ. The results show that the azimuthal 
deflections of myoV on actin are more common for the truncated lever arm, but the range of these deflections is 
relatively independent of its lever-arm length. 



INTRODUCTION 

Myosin V (myoV) is an unconventional myosin that ex- 
ists as a dimer in cells and is able to step processively 
along filaments of actin for relatively large distances 
(~1 jim) through a cytoplasm that is densely packed 
(Mehta et al., 1999; Reek-Peterson et al., 2000; Sakamoto 
et al., 2000; Walker et al., 2000; Veigel et al., 2002; 
Sweeney and Houdusse, 2010). How myoV is able to 
navigate through such a complex environment has 
been the subject of several in vitro studies that focused 
on dissecting the basic stepping mechanism of individ- 
ual myoV motors on isolated actin filaments (Rief et al., 
2000; Ali et al., 2002, 2007; Forkey et al., 2003; Toprak 
et al., 2006; Kodera et al., 2010). myoV, as with other 
members of the myosin family, binds in a fixed orientation 
to filamentous actin (F-actin) through its N-terminal 
motor domain. At the C-terminal end of the motor do- 
main is the light chain domain (LCD; also known as the 
lever arm), which is defined by six calmodulin (CaM)- 
binding IQ motifs, each bound to a CaM or CaM-like 
subunit. This domain amplifies small motions in the 
motor head associated with each ATP turnover result- 
ing in forward motion along the actin filament track, a 
concept known as the swinging lever-arm hypothesis 



Correspondence to Yale E. Goldman: goldmany@maiI.med.upenn.edu 

Abbreviations used in this paper: BR, bifunctional rhodamine; CaM, 
calmodulin; DOPI, defocused orientation and position imaging; EM, 
electron microscopy; F-actin; filamentous actin; LCD, light chain domain; 
myoV, myosin V; PFI, polarized fluorescence intensity; polTIRF, polarized 
total internal reflection microscopy. 



(Cooke, 1986). The LCD has been found to have addi- 
tional degrees of freedom, giving the myosin conforma- 
tional flexibility (Dobbie et al., 1998; Corrie et al., 1999; 
Veigel etal., 2002). 

Because F-actin is a two-stranded right-handed helix, 
there is a potential for myoV to take off-axis steps that 
would require a high level of flexibility. Electron micro- 
scope images of myoV frozen during processive motility 
along F-actin indicate that its two heads bind to actin 
predominantly spanning 13 actin subunits, with smaller 
subpopulations spanning 11 and 15 subunits (Walker 
et al., 2000). The distance spanned by 13 actin subunits 
corresponds to the 36-nm half-helical repeat of F-actin, 
consistent with the 36-nm step size of myoV (Mehta 
et al., 1999) . Two actin monomers 13 subunits apart have 
approximately the same azimuthal orientation (Amos 
and Amos, 1991 ) , leading myoV to walk relatively straight 
by stepping typically 13 subunits (^36 nm) per step. 

This view is supported by many single-molecule ex- 
periments in which myoV walks along actin bound to 
a rigid glass surface, which inhibits large azimuthal 
changes around the actin (Forkey et al., 2003; Toprak 
et al., 2006; Kodera et al., 2010). When myoV was ob- 
served to walk along filaments suspended away from the 
glass surface, though, the motors either walked straight 
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or with a gradual left-handed pitch, suggesting that the 
leading head of myoV lands on the 13th and, some- 
times, on the 11th subunit from its trailing head (Ali 
et al., 2002) . Twirling assays, which monitor the rotation 
of actin gliding along multiple motors fixed to a sur- 
face, also found that myoV walks straight along actin or 
it follows a gradual left-handed helical path (Beausang 
etal., 2008a). 

Polarized total internal reflection microscopy (pol- 
TIRF) and defocused orientation and position imaging 
(DOPI) measurements have also shown myoV walking 
relatively straight (Forkey et al., 2003), with some sud- 
den sideways tilts (Syed et al., 2006). polTIRF and DOPI 
are methods that directly measure the three-dimensional 
orientation of fluorescent rhodamine probes bifunc- 
tionally attached to a CaM (BR-CaM) bound to one 
of the myoV LCDs. In polTIRF, several polarized light 



intensities are detected and used to estimate the polar 
and azimuthal angles /3 F and a P , which describe the 
orientation of the probe's absorption and emission 
dipoles relative to the long axis of actin (Fig. 1 A) 
(Forkey et al., 2003). In DOPI, the shape of a slightly 
out of focus image provides the orientation of the 
probe (Toprak et al., 2006). These techniques have 
clearly demonstrated that the LCD switches between 
two distinct orientations associated with leading and 
trailing positions of the heads, providing strong sup- 
port for a "hand-over-hand" stepping model (Forkey 
et al., 2003). 

polTIRF and DOPI measurements report angles of the 
fluorescent probe. Because the rhodamine probe is not 
generally aligned with the myoV LCD, measurements of 
the probe orientation do not directly correspond to the 
orientation of the lever arm. In this paper, however, we 
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Figure 1 . Schematics of relevant frames 
of reference and polarized intensities. 

(A) Definitions of the orientation of the 
vector representing the rhodamine probe 
absorption and emission dipoles, /3 P , a P , in 
the Actin Frame (1; x A , Ja, za) and 0 P , 4>f, 
in the Lever Frame (2; x L , y L , zl) , which is 
itself shown, /3 L , ct h , in the Actin Frame (3). 

(B) Laboratory frame: zi ab „ atl>r j is aligned in 
the optical axis toward the microscope ob- 
jective; x, abmtory and jMoratory are parallel to 
the slide surface. Note that the zi^a/ory and 
ya axes, which are identical, point down- 
ward (Materials and methods) but are 
drawn pointing upward in A for clarity. Ex- 
perimental setup for polLIRF microscopy: 
A 532-nm laser is initially split into paths 1 
and 2 and is polarized using an optical setup 
similar to the one described in Beausang 
et al. (2008a). At the sample, the two beams 
in orthogonal planes, xiaboratory-^iaboratary and 
yiaimatarr^utimatmp alternatingly strike the 
quartz slide at angles (measured from 
the optical axis) slightly larger than the 
critical angle required for LIRF (Hecht, 
2001 ) . Each beam is sequentially polarized 
by Pockels cells into four polarizations: 
(1) horizontally (si and s2, perpendicular 

tO the ^laboratory ^-laboratory and yiaboratory ^-laboratory 

scattering planes for paths 1 and 2, re- 
spectively); (2) in the scattering planes 
(pi and p2, respectively); and (3 and 4) 
±45° (Ll/Rl and L2/R2, at ±45° angles 
with respect to the scattering plane). 12- 
and 16-channel measurements: Emission 
from the excited sample is split into its 
xiabaratory and yiaboratory components for detec- 
tion by avalanche photodiodes (Materials 
and methods), leading to a maximum of 
16 unique combinations of the eight time- 
multiplexed input paths and polarizations 
and two simultaneously detected polarized 
emission intensities. 12-channel measure- 
ments were made when the L and R polar- 
izations were not used in path 2. 



102 



Sideways motions of myosin V-6IQ and V-4IQ 



developed a way to determine the local angle of the 
probe relative to the lever arm and to use that informa- 
tion to calculate the actual lever-arm orientation. 

The step sizes of recombinant myoV constructs with 
different lever-arm lengths have been shown to depend 
direcdy on the length of their LCDs (Purcell et al., 2002; 
Sakamoto et al., 2005). In this work, the orientation of 
the BR-CaM-labeled lever arm of a truncated (4IQ) and 
full-length (6IQ) construct was compared using polTIRF 
microscopy. The truncated construct, indeed, walked 
much less straight along actin, wobbling from side to side 
as it progressed. We also determined the influence of the 
lever-arm length on the polar rotation ("twirling") of an 
actin filament translocated by multiple myosins in a glid- 
ing filament assay. The truncated construct led to more 
frequent actin twirling, surprisingly demonstrating either 
left- or right-handedness. 

MATERIALS AND METHODS 

Preparation of proteins 

G-actin was obtained from rabbit skeletal muscle (Pardee and 
Spudich, 1982). Biotinylated Alexa Fluor 647-labeled F-actin was 
prepared from G-actin, Alexa Fluor 647 actin (Invitrogen), and 
biotin-actin (Cytoskeleton) in a 5:1:1 ratio with 1 pM of total actin 
subunit concentration and stabilized with 1.1 pM phalloidin (In- 
vitrogen). 0.3% rhodamine-labeled F-actin was prepared from 
unlabeled actin and 6'-IATR rhodamine-actin (Corrie and Craik, 
1994) and stabilized with 1.1 pM phalloidin. Recombinant 
chicken myoV, with its full-length lever arm (myoV-6IQ; amino 
acids 1-1099) or a truncated lever arm (myoV-4IQ; amino acids 
1-863) and with a FLAG affinity tag at its C terminus, was co- 
expressed with CaM in SF9 cells (Fig. SI) (Purcell et al., 2002). 
Chicken CaM, with residues Pro66 and Ala73 mutated to cysteine, 
was expressed in Escherichia coli (Putkey et al., 1985) , purified, and 
labeled with bifunctional rhodamine (BR; provided byJ.E.T. Corrie, 
National Institute for Medical Research, Mill Hill London, UK; 
Corrie et al., 1998). myoV was labeled by exchanging endoge- 
nous CaM with exogenous BR-CaM at low stoichiometry (Forkey 
etal., 2003). 

Buffers 

M5 buffer (pH 7.6) contains 25 mM KC1, 20 mM HEPES, 2 mM 
MgCL, and 1 mM EGTA in deionized water. M5 + buffer is M5 plus 
10 mM dithiothreitol (DTT) and 100 pg/ml wild-type CaM (ex- 
pressed in bacteria; Putkey et al., 1985). The motility buffer for 
single-molecule myoV motility assays is M5 + buffer plus 4 pM ATP 



and 50 mM DTT. The motility buffer for actin twirling assays is 
M5 + buffer plus 100 pM ATP, 10 mM phosphocreatine (Sigma- 
Aldrich), 0.3 mg/ml creatine phosphokinase (prepared daily 
from powder; Sigma-Aldrich), and 50 mM DTT. 

Experimental apparatus 

In the original polTIRF setup developed in our laboratory, four 
time-multiplexed laser polarizations, each predominately aligned 
along the laboratory x, y, or z axis, illuminated the sample with 
two intersecting beam paths. The fluorescence emission was di- 
rected through a polarizing beam-splitting cube and onto two ava- 
lanche photodiodes (giving polarized detection along the x and y 
directions over an area corresponding to a 2-pm diameter disc at 
the sample) resulting in eight measured polarized fluorescence 
intensities (PFIs) that were used to estimate the probe orientation 
during every 40-ms excitation cycle. Symmetric orientations of 
the probe related by reflections across the Cartesian planes re- 
sulted in the same intensities, leading to an eightfold ambiguity 
in probe orientation (i.e., angles were limited to an octant of a 
sphere). In this work, additional laser polarizations aligned at 
±45° between the original polarizations were included (Beausang 
et al., 2008a), either in one beam path for a total of 12 PFIs or in 
both beam paths for a total of 16 PFIs, measured in 80-ms time 
intervals (Fig. 1 B) . The angular range determined in these two 
cases improved to a quarter-sphere and a hemisphere, respec- 
tively. Dipole symmetry of the probe prevents further increase in 
unambiguous angular detection range for polTIRF and polarized 
fluorescence methods in general. 

As described previously (Forkey et al., 2000, 2003, 2005; Beausang 
et al., 2008a), polarization data from single myosin molecules 
was recorded and used as input to an analytical model of the 
probe's fluorescence emission to determine its orientation (d laboratory, 
§ laboratory/ in the laboratory frame of reference (xi ab0 ratory, )u^, 
^■laboratory) , which is defined in Fig. 1 B. In general, orientations in 
this work are vectors represented by two polar coordinates; i.e., 
CI = (polar angle, azimuthal angle) , where the polar angle is defined 
from the +z axis to the vector, and the azimuthal angle is defined 
to the projection of the vector in the xy plane from the +x toward 
the +y axis. The probe and lever orientations in the various refer- 
ence frames are summarized in Table 1 . 

To compare molecules moving on actin filaments aligned in 
different directions in the sample plane, we transformed the 
orientation of the probe at each time point from the laboratory 
frame to a frame of reference based on the actin filament (x A , y A , z A ) 
(Fig. 1 A) (Forkey et al., 2000, 2003, 2005; Beausang et al., 
2008a,b). The z A axis is aligned with the barbed (plus) end of 
the actin filament, which is inferred from the motion of the 
myosin before recording the polarization data (Beausang et al., 
2008a) , y A is aligned with the optical axis pointing into the ob- 
jective, and x A lies in the sample plane perpendicular to z A and y A . 
Because of the dipole symmetry, the orientation of the probe 
in the actin frame is described equivalently by 0 P = (f} P , ct P ) and 



TABtE 1 

List of variables 



Variable 

(x A ,y A , z A ) 
flp = (Pp, <*p) 

fi'p= (180° - p P ,a P ± 180°) 
Al = (Pl, ol l ) 

(x L , ;)'£, Z L ) 

(0p, Vp) 

{ ■^laboratory) yiaboratory-i ^-laboratory) 
laboratory^ IP laboratory) 



Definition 

Axes describing actin frame of reference 
Probe orientation in the actin frame 
Probe orientation in the actin frame due to dipole symmetry 
Lever orientation in the actin frame 
Axes describing lever frame of reference 
Probe orientation in the lever frame 
Axes describing laboratory frame of reference 
Probe orientation in the laboratory frame 
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O'p = (180° — j8p, a P ± 180°). These two orientations completely 
describe the probe measured with 16 PFIs, but 12-PFI data re- 
sults in probe orientations that have an additional degeneracy 
resulting in four equivalent orientations. For convenience, the 
orientations of the 12 PFIs most consistent with the 16-PFI data 
are presented. Dwell periods for (fi P , a P ) are defined between 
two consecutive transitions in f} P . 

Single-molecule motility assay 

A precleaned fused silica slide (Corning HPFS grade; Quartz 
Scientific, Inc.) was further treated in an ion plasma cleaner for 
5 min and then spin coated with 2 mg/ml polymethyl methacry- 
late (PMMA; Sigma-Aldrich) in methylene chloride. The PMMA- 
coated slide was assembled into a 10-20-pl flow chamber with a 
glass coverslip and double-sided adhesive tape. Actin was adhered 
to the slide and aligned approximately with x^,^ by successive 
additions and 1 min of incubations of 1 mg/ ml biotinylated BSA 
(Sigma-Aldrich), 0.5 mg/ml streptavidin (Sigma-Aldrich), and 
100 nM biotinylated Alexa Flour 647-labeled F-actin, each fol- 
lowed by washes with M5 + buffer. myoV, containing BR-CaM, was 
introduced into the sample chamber at 10-1,000 pM in M5 + buf- 
fer, which contains 4 uM MgATP. 

Only molecules with clearly changing polarization recordings 
from a single probe, as indicated by step photobleaching to the 
fluorescence intensity baseline, were chosen for analysis, corre- 
sponding to 298 and 378 recordings for myoV-6IQand myoV-4IQ, 
respectively. After determining the maximum likelihood orienta- 
tion for each cycle of these recordings (described below) , only 
molecules with at least two clear orientational steps were retained 
for further analysis (n = 73 molecules [25%] for myoV-6IQ; n = 49 
[13%] for myoV-4IQ). The percentage of myoVs stepping twice 
before the photobleach is 1.9-fold higher for 6IQ than for 4IQ, 
which is likely related to the 1 .9-fold higher run length of myoV-6IQ 
reported earlier (Sakamoto et al., 2005). 

Actin twirling assay 

A flow chamber was assembled as above, except that the fused 
silica slide was coated with 2.5 mg/ml poly-L-lysine (Sigma- 
Aldrich) instead of PMMA, and then 0.2 mg/ml anti-FLAG anti- 
body (Sigma-Aldrich) was added and incubated for 1 min. 20 pi 
of M).2 mg/ml of unlabeled myoV-6IQor myoV-4IQ (i.e., without 
BR-CaM) was flowed in and incubated for 8 min. Exposed poly-L- 
lysine and fused silica were blocked with 2 x 20-pl washes of 
0.5 mg/ml BSA. Two 20-pl aliquots of pre-sheared, unlabeled F-actin 
(5 pM) were incubated for 1 min each to block any inactive myo- 
sin heads, and then unbound actin was flushed out with 20 pi of 
2 mM ATP in M5 + buffer. After two further washes of M5 + to remove 
residual ATP, 5 nM actin filaments, sparsely (0.3%) labeled with 
rhodamine, in M5 + buffer were incubated for 1 min, and then 
motility buffer was added to initiate filament gliding. Moving fluo- 
rescent spots were selected for polTIRF recording, and then the 
probe orientation was determined as described above. 

Theory and analysis for transforming probe angles 
to lever angles 

Model of myoV used in the analysis. myoV is composed of two 
heavy chains bound together through a coiled-coil region located 
past the C-terminal end of their lever arms (Reek-Peterson et al., 
2000). Our simplified model of myoV treats each half of the mole- 
cule as an actin-binding motor domain connected via a compliant 
joint to a rigid, inextensible lever arm, which is, in turn, connected 
to the other half of the molecule via a free swivel at the lever-tail 
junction (Fig. 1 A) . The orientation of the lever arm relative to the 
actin filament fl L = (j3 L , ct L ) is represented by a vector pointing 
away from the motor domain toward the C terminus. The orienta- 
tion of the probe relative to the lever-arm H P L = (0 P , q> P ) is defined 
in the lever (or local) frame of reference (x L , y L , z L ) (Fig. 1 A). 



The z L axis of the lever frame is aligned with the lever-arm vector; 
the x L axis is aligned along a vector defined by two residues (from 
C266 to F262) separated by one turn in an a helix located in the 
motor domain that is roughly perpendicular to the actin filament. 
The y L axis is then defined as perpendicular to both the x L and z L 
axes, and its direction is given by the cross product: y L = z L X x L . 
This definition assumes that the twist of the lever arm around its 
long axis (the y L angle; Corrie et al., 1999) is constant. Values for 
(6 P , <p P ) depend on the orientation of the probe bound to the 
CaM and on how the CaM is bound to the lever arm. Because the 
probe is rigidly bound to the CaM and the CaM is rigidly bound to 
the lever arm, we assume that the values for (8 P , <p P ) are constant 
for each molecule. 

A prediction of the hand-over-hand model is that the f} P and /3 L 
angles switch between two states as the lever arm transitions back 
and forth between its leading and trailing positions. To the extent 
that the two lever arms are rigid, f3 LiLeading + ji^Traiiing = 180°. This 
analysis assumes these two features. We do not necessarily expect 
(or observe; Fig. 4 B) only two angular values for a P and a L , as 
the azimuthal position of myoV depends on (a) the random ori- 
entation of the actin subunit to which myoV initially binds, and 
(b) any azimuthal changes myoV makes as it steps along the 
actin helix. 

The change in a P after two steps of myoV equals the corre- 
sponding change in a L ( 2 Aa P = 2 AaJ. Fig. 2 is a diagram of a la- 
beled molecule of myoV stepping straight along actin, showing that 
even when the myosin walks straight, the azimuthal orientation of 
the probe relative to the actin (a P ) changes after a step because it is 
not, in general, located within the plane of rotation. However, a P 
returns to its original orientation after a second step because the 
structure is identical, except for the translocation. In this simple 
(straight-walking) case, it is clear that the change in a P and a L over 
two steps (i.e., 2 \a P and 2 Aa L ) equals zero. For more compli- 
cated steps along actin, even with changes in a L , the two quan- 
tities are related by 2 Aa P = 2 Aa: L (Section 3 of the Appendix and 
Eq. A.37), if the change in j3 P over two steps ( 2 A/3 P ) is zero. In 
practice, | 2 Aj3p| < 10° is used. In summary, the assumptions of the 
analysis, 2 A/3 P = 0 and 2 \a P = 2 Aa L , are consequences of the 
hand-over-hand mechanism. f3 L ,uadm g + PL,Tmiiing= 180° applies if 
the lever arms are sttaight. 

Transforming the probe angles to lever angles. To directly trans- 
form the orientation of the probe in the actin frame ((3 P , ct P ) into its 
corresponding lever orientation (p L , a L ), the orientation of the 
probe relative to the lever-arm (0^ <p P ) has to be determined (Fig. 1 A 
and Sections 1 and 2 of the Appendix) . Through a series of rotation 
operations that align the actin frame to the probe frame (Fig. S2) , the 
probe's vector (v A ) and orientation (jSp, a P ) in the actin frame are 
given by (Section 1 of the Appendix): 



V ^ = ( sin a L ■ cos (p p + cos a L - cos P L ■ sin (p p ) ■ sin 6 p + cos a L ■ sin P L ■ cos 0 p 

V ^ = (— cos a L ■ cos (p p + sin a L - cos P L ■ sin (p p ) ■ sin 6 p + sin a L - sin P L • cos 6 P ( 1 ) 

V A z = -sin/? L ■ sin0 p ■ sin<p p + cos/? L ■ cos6 p 



and 

p P (e p ,<p p ,p L ) = K OS (v A J 

a p (B p ,<p p ,a L ,P L ) = atan2(w t- , t ^ j ) 

The atan2(x, y) function is similar to atan (ji/jc) , except that atan2 
(x, y) is single valued over a larger range (-» ^ a F < n). 
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Inverting Eqs. 1 and 2 gives the lever-arm vector (w A ) and orien- 
tation (fl L , a L ) in the actin frame (Section 2 of the Appendix): 



with the following definition for the weighting function, W (Sec- 
tion 4 of the Appendix) : 



w A,x =(c°s(PA e r->P r -Pr))- sin <P r sinfl , +sin(j3 t (e p ,?> F ,j8,))-co s 8,) 
•cosa p — sin9 p ■ cosq> p ■ sina p 

W At y=(co S (pJe p ,<p r ,P r ))sin q , r sine p + S m(P L (0 p , Vr ,l3 p ))-cose r ) 
■sina p + sinf^ ■ cos(p p - cosa p 

W A,z = ( cos # P ±^jsm0 p 2 -sinfp/ -(sin/?/ -sm0 p 2 -cos<p/)j 

/ cosO p z + sinG p 2 -sintpp 2 



IV', 



and 



a L {e p ,<p p ,a p ,P p ) = atan2 (w Ax , w Ay ) 



(3) 



(4) 



Eqs. 3 and 4 imply that calculating j3 L and a L from the measured 
values of /3 P and a P requires finding reliable values of (8 P , <p P ) for 
each observed myoV. 

Determining 0 P and <p P . Differences in probe and lever-arm azi- 
muth after one step are given by Act p „ = a p „ + j — a Pn and Aa Ln = 
a L,n+i ~ a L.n< respectively, where n indexes the steps (n= 1, 2...iV— 1) 
between the N pauses in a processive run (Sections 3-5 in the Ap- 
pendix) . These two expressions are related by: 



{\ Act F,„ 



(10) 



Aa„ + e 



where a small value, e, is chosen to avoid division by 0, and ris an 
even positive integer. The results of the procedure are not sensitive 
to the values of e and r. Setting them to s = 0.0001 rad and r= 4 
results in a robust stable value of Rior each molecule. 

Given JV — 1 values of Aa P measured from a molecule taking N 
steps along actin, Eq. 9 is used to calculate an approximation for R 
This value for R is then applied in Eq. 6 to calculate the N— 1 cor- 
responding values for Aa L . 8 P and <p P are then adjusted in Eqs. 7 
and 8 using the N— 1 values of Aa L and T °'& L ,n = Pl,„ + j3 L ,„ +1 = 180°, 
assuming that the lever arm is straight in both the leading and 
trailing positions. Once an estimate for (flp, <p P ) is obtained for that 
molecule, /3l and a L are calculated using Eqs. 3 and 4 and the 
AT measured values of /3 P and a P . 

Treatment of degeneracy of probe orientation for myoV under- 
going multiple p P transitions. When using 16 channels of polar- 
ized intensity, the solution of the probe's orientation 0 P = (/3 P , a P ) 
is limited to a hemisphere of orientation, which is set to be 0° < p P < 
90° and — 180° < a P < 180°. Because of the symmetry of the dipole 



l ^P,n = A «P,„ + A «P,„ + 1 = A «L,„ + A «L. n+ l = 2A «L.„. (5) 



and thus 



Aa„ 



Aa, 



■ Aa, 



■Aa t 



= (-l) ml (Aa„ ,, -Aa, i( ,,) = R, 

V / \ P,n+{m-L) L,n+{m-L) J ' 



(6) 



where m is an integer describing the interval of steps being consid- 
ered, and R is a constant for a given molecule that depends on the 
local probe angle (0 P , <p P ) (Section 4 of the Appendix). When the 
molecule walks straight (Aa L = 0), R= |Aa: P \. 

Our strategy for calculating (9 P , <p P ) requires finding an estimate 
for R so as to calculate the N — 1 values of Aa^ n to use in the fol- 
lowing system of equations based on Eqs. 3 and 4 (Section 5 of the 
Appendix) : 

a L (e p ,<p p , a p iB+ , , p p B+1 )-a L (0 p ,<p p ,a P n , p pn ) = Aa L „ 



PA°»<P» ft™ ) + ft ft,, ) = m p L ,„ ■ 



(») 



The value for R for each molecule is obtained by calculating the 
average of (— l)"-Aa: P over all n. If ka L is zero or a constant for a 
processive run, this average exactly equals R (Section 4 of the Ap- 
pendix). \a L (and by extension, | Aa: P |) is not generally constant, 
so a weighted average of ( — l)"-Aa: P (i.e., Aa p ) is calculated that 
deemphasizes values of Aa P that occur when the molecule takes 
large sideways steps that are atypical: 
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Figure 2. Schematic explanation of 2 Aa. The azimuthal changes 
in probe angle, a P , and lever angle, a L , over two steps are referred 
to as 2 Aa: P and Aat, respectively. When walking straight, the probe 
attached to the myoV lever arm returns to its original orientation 
every two steps, so that both 2 Aa P and 2 Aa L equal zero and 2 Aa P = 
2 Aa t . It can be shown that 2 Ao: P = 2 ha L is a general relationship, 
regardless of straight walking, so long as f3 P switches regularly be- 
tween two states (Section 3 of the Appendix). 
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(fl'p = (180° — p P , a P ± 180°)), the probe has an equivalent orien- 
tation in the hemisphere defined by 90° < fi P < 180° and —180° < 
a P < 180° (Forkey et al., 2000). As mentioned in the last section, 
Pp regularly switches between two angular states as myoV takes 
N — 1 steps along actin. These two states, along with the twofold 
dipole symmetry, lead to four equivalent expressions for the ^mea- 
sured orientations of the probe. We express these four equivalent 
orientations as N — 1 pairs consisting of two sequential angular 
states: (n Pn , f2 Pn+ i), (f2 Pm fl'pn+i)) (ft p jIU f2p,n+i)> and (ftpm ftp, n +i)' 
These four equivalent probe pairs can be used to calculate four 
corresponding pairs of the lever-arm orientation for the N angu- 
lar states of n Ln = (fi^ a Li „) = (P L (n P , m 9 P , <p P ), ot L (n P , n , 6 P , <p P )) 
(Eqs. 3 and 4) . Solutions for a L „ are the same for these four angu- 
lar paths, but for /3 i n , there are two solutions that result in the 
different lever-arm configurations, A and B (Fig. S3). Lever-arm 
configuration A is calculated from the probe angle pairs (ftp,„, 
fip n+ i) and (0'p n , fi'p n+ i) using the local probe orientations (B P , 
<p P ) and (180° — Bp, <p P ± 180°), respectively (Fig. S3 A); configura- 
tion B is similarly calculated from the probe angle pairs (f2 Pn , 
fi'pn+i) and (fi'p n , n PiI1+1 ) (Fig. S3 B). In general, there is not a 
straightforward relationship between these two lever-arm configu- 
rations except for the relationship between the local probe angles 
(fip and O'p). 

For each of the four configurations described above, there is an 
additional equivalent solution for the lever-arm orientation that 
corresponds to their mirror images (Fig. S4). Given any probe 
solution Hp;, two equivalent lever solutions (n L , n = flY.n) can 
be calculated that satisfy Eqs. 7 and 8, where 0 Ln is described 
above and 0' L ,n = (180° - /3 Ln , a hn ± 180°) = (p L (ftp,n, 0' P , <p' P ), 
a L (n P?n , 0'p, <p' p )), where 6' P = 180° - 9 P , <p' P = 180° - (p P . Note that 
this symmetry is independent of the dipole symmetry of the 
probe; so with this mirror symmetry and the four lever solutions 
related to the dipole symmetry, there are a total of eight possible 
related solutions for the lever-arm orientation. 



Online supplemental material 

Fig. SI provides a basic schematic of the myoV-6IQ and myoV-4IQ 
recombinant constructs. Fig. S2 describes the different rotation op- 
erations performed by the rotation matrix used in deriving Eqs. 1 
and 2 (Section 1 of Appendix). Figs. S3 and S4 demonstrate the 
degeneracy of the lever angle solution (Materials and methods). 
Figs. S5, S8, and S9 give various probe and lever angle distributions 
for myoV-6IQ measured using 16 PFIs (Materials and methods). 
Tables S1-S3 give the best-fit values from the Gaussian fits to these 
distributions. Figs. S6 and S7 give the distributions for Aa P and 2 A/8 P 
for myoV-6IQ and myoV-4IQ. Table S4 gives the values from the 
Gaussian fits to these distributions. Fig. 10 illustrates the twirling 
experiment. Fig. SI 1 and Table S5 compare the polTIRF data from 
myoV-6IQ and myoV-4IQ to a molecular model of myoV (Parker 
et al., 2009) . The supplemental material is available at http:/ /www 
.jgp.org/cgi/content/full/jgp.201110715/DCl. 

RESULTS 

Measurements of the three-dimensional orientation 
of myoV-6IQ and myoV-4IQ during single-molecule 
processive motility 

polTIRF measurements were made with BR probes 
(Forkey et al., 2003, 2005) bound to CaM subunits on 
two recombinant constructs of myoV derived from 
chicken: myoV-6IQwith its native number of CaMs (six) 
bound, and myoV-4IQ, truncated to contain four CaMs 
(Purcell et al., 2002). Both constructs have a coiled-coil 
region at the C terminus for dimerization but no cargo- 
binding domain (Fig. SI ) . Probe angular measurements 
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Figure 3. Representative traces from the 
processive runs of rhodamine-labeled myoV- 
6IQ (left) and myoV-4IQ (right) motors 
along actin filaments. (A and B) Measured 
intensities: The total measured intensity 
is shown (A), along with its 12 compo- 
nent intensities (B). When the rhodamine 
probe bleached, all of the intensities 
dropped to baseline. (C-E) Calculation of 
the orientation of the probe in the actin 
frame of reference: Traces of the angles 8 P , 
microsecond wobble (C), f} P , probe axial 
angle (D), and a P , probe azimuthal angle 
(E) in the actin coordinate frame. 
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for most of the data were made using 12 PFI channels 
(Figs. 1 B and 3-8, and Tables 2-4); further data for 
myoV-6IQwas taken using 16 PFI channels (Figs. S5-S9 
and Tables S1-S4) . Full sets of these PFIs were recorded 
at 80-ms intervals. 

Video sequences of BR on single myoV molecules 
moving processively on actin attached to the micro- 
scope slide gave average translocation velocities, V ava of 
103 ± 3 nm/s (mean ± SEM; n = 92) and 85 ± 3 nm/s (n = 
50) for myoV-6IQ and myoV-4IQ, respectively, at 4 uM 
MgATP, consistent with prior findings (Cheney et al., 
1993). Representative recordings of the 12 individual 
polTIRF intensities (Fig. 3) for myoV-6IQ, and myoV-4IQ 
stepping at 4 pM MgATP show that the individual traces 
are typically stable for periods of several hundred 
milliseconds and then suddenly transition, back and 
forth, to new intensity levels (Fig. 3 B). Complementary 
changes among the 12 traces result in the total intensity 
remaining relatively constant (Fig. 3 A) . The dwell peri- 
ods were exponentially distributed with mean dwell 
times of 702 ± 23.7 ms (mean ± SEM) and 927 ± 53.6 ms 
for myoV-6IQ and myoV-4IQ, respectively. The polar 
(Fig. 3 D, (3p) and azimuthal (Fig. 3 E, a P ) angles of the 
BR probe relative to the actin, calculated from the pol- 
TIRF data as described in Materials and methods, show 
regular switching between two orientations. The wobble 
angle, 8 P , is noisy (Fig. 3 C) but relatively stable. These 
features of the polTIRF recordings are consistent with 
hand-over-hand motility generated by tilting lever arms 
in both constructs. 

Probe orientation distributions for myoV-6IQ 
and myoV-4IQ are similar 

The traces of the axial and azimuthal angles of the 
probe, /3p and a P , respectively, are stable during dwell 
periods. Mean values of /3p and a P , calculated over these 
stable intervals, are plotted as histograms, where the (3 P 
distribution is fit to two Gaussian components (Fig. 4) . 
The two components of the /3 P histogram, correspond- 
ing to the leading and trailing head positions, peak at 
80° and 25° for myoV-6IQ, respectively, and 81° and 31° 
for myoV-4IQ (Fig. 4 A and Table 2) . Note that an orien- 
tational ambiguity as a result of the symmetry of the op- 
tical dipoles limits expression of these probe angles to 
the range of 0°-90°, a restriction that does not apply to 
calculations of lever-arm angles given later. The values 
for /3p are similar to those reported previously (Forkey 
et al., 2003; Toprak et al., 2006) , when this orientational 
ambiguity is considered (Materials and methods) . 

a P usually visits two values during individual runs of 
myoV-6IQand multiple values for myoV-4IQ (Fig. 3 E). 
Distributions of a^from many processive runs of both 
constructs are broad (Fig. 4 B) , most likely because the 
starting azimuth depends randomly on which actin sub- 
unit the myoV initially binds. To better understand 
the nature of the azimuthal paths and to eliminate the 



variability generated by the starting subunit, we instead 
consider changes in a^as myoV steps along actin. These 
changes are independent of the starting azimuth. 

2 \a P and 2 A/3 P distributions for myoV-6IQ and myoV-4IQ 

The change of the azimuthal angle of the probe after 
one step of myoV, a Pn+1 — a Pn , is termed Aa^. (Fig. 2). 
One should note that even if the molecule walks straight 
along actin, Aa^will generally be nonzero because the 
BR probe is not located direcdy in the plane of lever- 
arm rotation. However, the quantity 2 Aa P = a Pr , +2 — a Pn , 
the change in a P after two steps, does reproduce the 
azimuthal change of the lever arm because in hand- 
over-hand motion, the molecule returns to the same 
state every other step, except for the axial translation 
and any accumulated azimuthal motion. 

The azimuthal angle of the probe, a B is related to that 
of the lever arm of myoV, a L , through the expression 
2 Aa P = 2 Aa L (Section 3 of the Appendix). This equality 
between ^Ao^and 2 Aa L is a general one, independent of 
the local orientation of the probe with respect to the 
lever arm (given by polar angle 6 P and local azimuth cp P ; 
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Figure 4. f} P and a P distributions for myoV-6IQ and myoV-4IQ. 
Histograms of the (3p and a P angles were collected from the mean 
angles in states that were visited during individual runs of myoV- 
6IQ (top: n = 73 molecules and states = 515) and myoV-4IQ (bot- 
tom: n = 49 molecules and states = 246) . The beginning and end 
of each dwell were defined by the transitions in /3 P . Histograms 
are normalized so that the total area of the distribution is equal 
to one. (A) f} P distributions: Each of these distributions was best 
fit by a bimodal Gaussian equation; the resulting parameters are 
given in Table 2. (B) a P distributions. 
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TABLE 2 

Analysis of the distributions of probe a and f3 angles for myoV-6IQ and 
myoV-4IQ in the actin frame of reference 



Variable 



3 Trailing 



2 Aa, 
2 Aa 2 
2 Aa 3 



myoV-6IQ myoV-4IQ 
Pp angles (°) 
80 ±12 (52%) 81 ±11 (51%) 

25 ± 8.0 (48%) 31 ± 8.3 (49%) 

2 Act P angles (°) 
-30 ±27 (31%) NA 
-2.0 ± 14 (50%) -9.7 ± 88 (100%) 

34 ±27 (19%) NA 



Values result from fitting a Gaussian model to the distributions and report 
peak ± (j. 

Fig. 1 A), with the single condition that the change in (3 P 
over two steps ( 2 A/3 P ) must be equal or near to zero, as 
expected in a hand-over-hand mechanism (Materials and 
methods). The measured 2 A/3 P distributions for myoV- 
6IQ and myoV-4IQ are narrow and centered on zero 
(Fig. S7 and Table S4) , indicating that this condition is 
satisfied in our measurements. The ~Aa P distribution for 
myoV-6IQhas one principle peak at —2.0 ±14° (peak ± a) 
corresponding to straight walking by the motor 13 
actin monomers per step. Two shoulders on either side 
of the central peak can be fitted by further Gaussian 
peaks at -30 ± 27° and +34 + 27° (Fig. 5 A and Table 2) , 
where the widths of the side components were con- 
strained to be equal. This three-peak distribution fits the 
data better than a single Gaussian component (P < 0.0001 , 
as determined by an F-test). The two side lobes in the 
myoV-6IQ "Aa P distribution correspond approximately 
to the motor stepping 11 or 15 subunits for one of the 
steps. Two Gaussian peaks (at -2.7 + 12° and -6.9 ± 39°) 
also fit the data. The central peak is the same and the 
overlapping wider component corresponds to the two 
side lobes discussed above. In contrast, the 2 Aa P distribu- 
tion of myoV-4IQ (Fig. 5 B) is much broader without 
distinct subcomponents. A single Gaussian peak, —9.7 ± 
88° (Table 2), describes the histogram of myoV-4IQ 2 Ao'/). 

Multiple subpopulations in the distributions of 0 P versus <p P 
for myoV-6IQ and myoV-4IQ 

The orientation of the BR probe relative to a frame of 
reference fixed within the lever arm is described by the 
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polar angles (6 P , cp P ) (Materials and methods and 
Fig. 1 A) . Using analytical relationships between the ori- 
entation of the probe (f3 P , a P ) and lever-arm {(5 L , a L ) 
relative to actin and (6 P , <p P ) (Eqs. 7 and 8), we can esti- 
mate (6 P , (p P ) from the jS^and a P angles measured from 
a series of steps of myoV in a processive run. This esti- 
mation arises from the relationship 2 Act P = 2 Aa L (when 
2 Af3 P = 0) and the further assumption that the two lever 
arms are rigid, leading to the relationship 7oi /3 L = ^L,TraiUng + 



2 Aa„ 



ml 



L,Leading ' 



180°, where f3 



L, Trailing 



and /8 



L,Leading 



are the axial 



angles corresponding to the two head positions. 

Scatter plots of flpversus ^calculated from individual 
runs of myoV-6IQ and myoV-4IQ each show three fairly 
distinct populations (Fig. 6 A) . Each of these groups was 
separately binned into two-dimensional histograms and 
fit by two-dimensional Gaussian functions, giving peak 
values of (0 P , cp P ) = (45°, -110°), (43°, -73°), and (63°, 
— 65°) for myoV-6IQ (Table 4). The distribution of 
calculated values of (6 P , <p P ) for myoV-4IQ was quite 
similar, also demonstrating three apparent populations 
peaking at (43°, -113°), (46°, -60°), and (65°, -62°) 
(Table 4). We tentatively hypothesize that these compo- 
nents of the (6 P , cp P ) distributions for the two constructs 
correspond to different CaM subunits that exchanged 
with BR-CaM during labeling. 

fi L and a L were then calculated from the measured 
values of (3 P and a P , the estimated (6 P , (p P ) for each mol- 
ecule, and Eqs. 3 and 4. Representative traces of f3 p and a L 
(Fig. 6 B) indicate that myoV-6IQ walks hand-over-hand 
along actin with little azimuthal deviation (a L is con- 
stant) , sporadically exhibiting sideways motions, as 
described previously (Syed et al., 2006). myoV-4IQ dis- 
plays similar values for (5 L to those of myoV-6IQ but has 
a variable a L , wandering azimuthally around actin as it 



walks. Distributions of j3 



L, Trailing 



and £ 



L,Leading 



values for 



myoV-6IQ ( n = 73 traces) and myoV-4IQ ( n = 49 traces; 
Fig. 7 A) show similar pairs of peaks for the two 
constructs. Distributions of a L for both constructs are 
broad, as expected again because of their random start- 
ing azimuthal positions (Fig. 7 B). 

There are two possible unique solutions for the lever 
angles that describe different lever-arm configurations, 
A and B (Materials and methods and Fig. S3). The dis- 
tributions of j5 L calculated for these two configurations 
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Figure 5. 2 JAa P distributions for myoV-6IQ and myoV-4IQ. 
Histograms of the 2 tAa P angles were calculated from the 
mean angles of a P that were visited during runs of myoV- 
6IQ (A) and myoV-4IQ (B) and normalized as in Fig. 4. 
The distributions were fit by a trimodal (6IQ) or a uni- 
modal (4IQ) Gaussian equation; the resulting parameter 
values are given in Table 2. 
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for myoV-6IQ (measured with 16 channels) each have 
two distinct peaks that are not substantially different 
(Fig. S5 A and Table S2). Therefore, although there 
are four equivalent probe orientations from a given 
run of a labeled myoV, each of these results in roughly 
similar values for f3i, thus reducing the final degener- 
acy of lever-arm orientation from eightfold to twofold. 
Of these two possible solutions, we eliminate one by 
assuming that there is not enough clearance for myoV 
to walk between the actin filament and the slide. Based 
on this assumption, we then accept the lever orienta- 
tions that place the molecule on the solution side of 
the actin filament away from the glass (i.e., 0° < ct L < 
180°) (Fig. S5B). 

Most of our data were measured using 12 channels of 
polarized intensity, which limits the probe to a quarter- 
sphere of space and leads to further degeneracies in the 
probe and lever-arm orientations, as mentioned earlier. 
We present a single solution for the 12-channel data for 
each molecule by constraining the probe and lever ori- 
entations to be close to those found for the 16-channel 
myoV-6IQ data. 

Aa L distributions for myoV-6IQ and myoV-4IQ 

Distributions of Aa L (= a Li , l+1 — a L n ), calculated from 
(a F , /3p) and (Op, (pp) for each molecule (Fig. 8), report 
the azimuthal changes of the lever arms after single 
steps of myoV along actin. The sign for Aa L is negative 



when myoV walks leftward, at or near zero when myoV 
walks straight, and positive when walking rightward. 
The distribution of Aa L values from myoV-6IQ data 
has a single principle peak at —1.8 + 13° and two ap- 
parent shoulders at -29 ± 34° and +33 ± 34° (Fig. 8 A 
and Table 3) , where the widths of the two side compo- 
nents were constrained to be equal. The straight-walking 
component of steps (Aa L = —1.8°) makes up ^55% 
of the total population. The Aa L distribution for 
myoV-4IQ_ also has three populations with peak values 
very similar to those from the myoV-6IQ_ distributions 
(Aa L = -35 ± 29°, -1.4 ± 10°, and +47 ± 29°; Fig. 8 B 
and Table 3) , where the widths of the side components 
were constrained to be equal. These three-peak dis- 
tributions fit the data better than single Gaussian 
components (P < 0.0001 for both 6IQ and 4IQ_, as 
determined by an F-test). As with the myoV-6IQ_ 2 Aa P 
distribution described earlier, two Gaussian-centered 
peaks with different widths (at -2.3 ± 46° and —1.9 ± 
12° [6IQJ; -3.3 ± 54° and -1.9 ± 8.7° [4IQJ) also fit 
the data. The central peaks are the same, and the over- 
lapping wider components correspond to the two 
side lobes. The straight-walking population from the 
myoV-4IQ_ distribution (Aa L = —0.2°) makes up only 
36% of the total population, much less than that for 
myoV-6IQ. Thus, an individual molecule of myoV-4IQ 
is more likely to take steps to the left or right than to 
walk straight. 
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Figure 6. 6 P versus (p P distributions for myoV-6IQ and 
myoV-4IQ. (A) 0 P versus <p P distributions calculated from 
individual processive runs of myoV-6IQ (left; n = 73 mol- 
ecules) and myoV-4IQ (right; n = 49 molecules) describing 
the orientation of the probe in the lever frame of reference. 
One-dimensional histograms of 8 P and <p P are displayed to 
the left and below them. (B) Representative traces of p L 
and a L , the orientation of the lever arm in the actin refer- 
ence frame, calculated from the /Jpand ct P traces shown in 
Fig. 3 along with the corresponding values for 6 P and <p P in 
those molecules (6IQ: 8 P = 62° and <p P = -65°; 4IQ: 0 P = 42° 
and (p P = -124°). 
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TABLE 3 

Analysis of the distributions of lever-arm a and fj angles for myoV-6IQ 
and myoV-4IQ in the actin frame of reference 
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Figure 7. fi L and a L distributions for myoV-6IQ and myoV-4IQ. 
Histograms of the /3 L and a L , mean angles visited during indi- 
vidual runs of myoV-6IQ (top: n = 73 molecules and states = 515) 
and myoV-4IQ (bottom: n = 49 molecules and states = 246). 
Dwell period selection and normalization as in Fig. 3. (A) f} L 
distributions: Each of these distributions was fit by a bimodal 
Gaussian equation; the resulting parameters are given in Table 3. 
(B) a L distributions. 



Twirling of actin in gliding assays with myoV-6IQ 
and myoV-4IQ 

Unlabeled myoV-6IQ_ or myoV-4IQ_ was attached to the 
microscope slide using a C-terminal Flag tag (Fig. SI) 
and anti-Flag antibody. Actin, very sparsely labeled with 
(monofunctional) rhodamine (Materials and meth- 
ods), glided on the myosin (Fig. S10) with velocity of 
112 + 5.8 and 133 ± 7.0 nm/s for myoV-6IQ and myoV- 
4IQ, respectively, at 0.1 mM MgATP. The angle of the 
fluorescent probe in actin was determined during glid- 
ing by recording 16 PFIs (Materials and methods) 



Variable 

Leading 
^Trailing 

Aa, (L) 
L\a 2 (S) 
L\a 3 (R) 



myoV-6IQ myoV-4IQ 

j8 £ angles (°) 
120 ± 12 (50%) 124 ± 15 (52%) 

60 ±13 (50%) 56 ±16 (48%) 

Acti angles ( 0 ) 
-29 ± 34 (25%) -35 ± 29 (36%) 

-1.8 ±13 (55%) -1.4±10(36%) 
33 ± 34 (20%) 47 ± 29 (28%) 



Values result from fitting a Gaussian model to the distributions and report 
peak ± cr. 



(Beausang et al., 2008a). As with the other myosins 
measured in the twirling assay, for both myoV-6IQ and 
myoV-4IQ, the polTIRF intensities gradually rose and 
fell, and the individual polarizations exhibited various 
relative phase shifts (Fig. 9 B) . Total intensity was con- 
stant until the probe photobleached, indicating that the 
changes of the individual polTIRF traces are caused by 
probe rotation as the actin twirled around its axis dur- 
ing gliding (Beausang et al., 2008a). /3p often changed 
very little over the course of the experiment, as expected 
from a probe fixed to the actin (Fig. 9 D) . On the other 
hand, a P decreased at a constant rate with myoV-6IQ_ 
and either steadily increased or decreased with myoV-4IQ_ 
(Fig. 9 E). For each actin filament, we calculated that 
the pitch = y/w of twirling, where V is the average 
linear velocity (nm ■ s _1 ) of gliding and oo is the angular 
velocity (rotations ■ s _1 ) calculated from Aa P /At. Inas- 
much as some of the filaments did not twirl (as was also 
found before for myosins II, V, and VI; Sun et al., 2007; 
Beausang et al., 2008a) , we calculated the inverse pitch 
(pitch 1 (um~ )) to display distributions of both twirl- 
ing and non-twirling filaments on single distributions 
(Fig. 10). Larger values for the inverse pitch indicate 
tighter twirling, and positive and negative signs indicate 
right- and left-handed twirling, respectively. 

myoV-6IQ and myoV-4IQ show two clear differences 
in their twirling. A significant portion of filaments 
gliding over myoV-6IQ-coated surfaces did not twirl 
(IpitcrT 1 ! < 0.4 pm = 50%), whereas most filaments 
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Figure 8. \a L distributions for myoV-6IQ and myoV-4IQ. 
Histograms of the Ao: £ angles calculated from the mean 
angles of a L that were visited during individual runs of 
myoV-6IQ (A) and myoV-4IQ (B). The histograms are 
normalized to one. Each of these distributions was fit by a 
trimodal Gaussian equation; the resulting parameters are 
given in Table 3. 
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TABLE 4 

Analysis of the 6 P versus (p P distributions calculated from the processive 
runs of myoV-6IQ and myoV-4IQ_ 



Probe position 


M°) 


<Pr (°) 


% Total 






myoV-6IQ 




1 


45.2 + 7.67 


-110 ±5.12 


38 


2 


43.2 ± 7.96 


-73.0 ± 12.8 


28 


3 


63.0 ± 4.05 


-64.9 ± 7.06 


34 






myoV-4IQ 




1 


42.9 ± 5.39 


-113 ±11.1 


26 


2 


45.9 ± 5.71 


-59.8 ± 8.83 


45 


3 


64.5 ± 6.59 


-61.9 ±4.45 


2!) 



Values result from fitting a Gaussian model to the distributions and report 
peak ± (T. 



on myoV-4IQ slides twirled ( |pitch 1 1 < 0.4 pm = 
6%) . From those filaments that did twirl, a left-handed 
pitch of 1.4 ± 0.13 um (mean ± SEM) was calculated, 
in agreement with earlier measurements on native 
myoV (Beausang et al., 2008a). When the pitches of 
both the twirling and non-twirling filaments from 
myoV-6IQ are considered, an average left-handed 
pitch of 2.7 ± 0.64 um is calculated. myoV-4IQ, sur- 
prisingly, twirls either to the left (56%) or to the right 
(38%), with pitches of 1.2 ± 0.14 pm and 1.0 ± 0.19 
pm, respectively. myoV-4IQ is the first myosin tested 
that can twirl robusdy with either handedness, suggesting 



a coupling between the step sizes and azimuthal changes 
among motors mechanically connected through actin 
and the glass slide. 

DISCUSSION 

Truncating the lever arm of myoV increases off-axis steps 
In this study, we were able to calculate the changes in 
orientation of the lever arm of myoV as it steps along 
actin using a novel analysis method. Previous single- 
molecule studies using polTIRF characterized the an- 
gular changes of fluorescent probes on the lever arm, 
denoted a P and j3p, associated with stepping by myoV 
(Forkey et al., 2003; Syed et al., 2006; Toprak et al., 
2006) , but they did not transform these angles into the 
actual orientation of the lever arm, a L and )3l, because 
the local orientation of the probe relative to the lever 
was generally unknown. The gradual azimuthal change 
as myoV travels over several steps (>1 pm) , as measured 
by its helical pitch, has been measured in other single- 
molecule (Ali et al., 2002) and ensemble studies 
(Beausang et al., 2008a), but these earlier studies did 
not quantify the step-by-step changes in its azimuth, as 
measured by Aa L . We were able to calculate distribu- 
tions for Aa L from the present measurements of the ori- 
entation of the probe attached to individual myoVs. 
The native myoV construct, myoV-6IQ, predominandy 
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Figure 9. Representative traces from 
rhodamine-labeled actin filaments twirling 
over myoV-6IQ (left) or myoV-4IQ (right). 
Measured intensities (A and B) : The total 
measured intensity is shown (A), along 
with its 16 component intensities (B). The 
rhodamine probe bleached when the in- 
tensities dropped to baseline. (C-E) Cal- 
culation of the orientation of the probe 
in the actin frame of reference: Traces of 
the angles S P (C), j3 P (D), and a P (E) in 
the actin frame. Note that a P declines lin- 
early for myoV-6IQ, indicating left-handed 
twirling, and that a P increases linearly for 
myoV-4IQ, indicating right-handed twirl- 
ing in this particular molecule. 
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takes straight steps, as indicated by the near-zero primary 
peak that makes up 55% of the Aa L distribution (Fig. 8 
and Table 3) . This high likelihood for taking a straight 
step agrees well with the results from the ensemble-level 
gliding assays, where 50% of the actin filaments did 
not twirl (i.e., they moved straight) (Fig. 10), as was ob- 
served previously for tissue-purified myoV (Beausang 
et al., 2008a). Small shoulders to the left and right of 
the principle Aa L peak, at -29° and +33° (Fig. 8), ac- 
count for the remainder of the Aa L distribution, where 
leftward steps (Aa L < 0) are found almost twice as often 
as rightward steps ( Aa L > 0) (Table 3) . This leftward 
stepping is in agreement with the behavior of the other 
50% of the observed filaments in the gliding assay, 
which twirled with a left-handed pitch (Figs. 9 and 10), 
as observed previously (Ali et al., 2002; Beausang et al., 
2008a). From the entire Aa L distribution, the average 
change in a L is —2.3° per step. With an average step size 
of 36 nm, the helical path is predicted to be left-handed 
with a pitch of 5.6 pm (360° ■ 36nm/2.3°). This value is 
approximately twofold larger than the overall average 
pitch we found in the ensemble gliding assays and what 
was reported by Ali et al. (2002) using an assay with sus- 
pended actin. This difference is likely to be caused by 
the geometry of the polTIRF system, in which actin is 
fixed to the surface, preventing myoVfrom stepping all 
the way around the actin filament and thereby altering 
the number of sideways steps. 

myoV-4IQ, with a truncated lever arm, has a Aa L dis- 
tribution with a central peak near zero and two shoul- 
ders to the left and right, all with means similar to those 
for the native myoV-6IQ (Fig. 8 and Table 3) . However, 
the two constructs differ significantly in the relative fre- 
quencies of these three subpopulations; myoV-4IQ has 



an approximately equal probability for stepping left, 
right, or straight for any given single step, whereas 
myoV-6IQ walks straight more often. These results are 
consistent with those of the myoV-4IQ gliding filament 
assays, where most of the observed actin filaments 
twirled (94%) with a left-handed (56%) or right-handed 
pitch (38%), a loss of azimuthal bias in the twirling of 
actin that has not been observed before for any other 
myosin family member (Sun et al., 2007; Beausang et al., 
2008a). In the kinesin family of microtubule motors, 
a somewhat analogous observation has been made 
for Ned, where a point mutation in its neck linker 
region led to a loss of directional bias in microtubule 
gliding assays (Endow and Higuchi, 2000). Individual 
microtubules were driven by this mutant in either the 
plus- or minus-end direction. Such behavior in the 
ensemble assay suggests that cooperation between the 
motors that are linked through the cytoskeletal filament 
promotes directional concordance of individual motors 
that have limited directional bias, thereby maintaining 
an initial directionality. In the case of myoV-4IQ ob- 
served here, the axial direction along the actin filament 
(toward the barbed end) is not variable, but the fre- 
quency for an individual motor to step sideways is high, 
as discussed below, and leftward and rightward steps 
have similar probabilities. Linking the motors through 
the filament apparently coordinates this helical direc- 
tion, leading to either left- or right-handed twirling. 

Correlating the azimuth of stepping myoVto the helical 
structure of F-actin 

Subunits of F-actin form a right-handed helix with a he- 
lical half-pitch repeat of ^13 subunits that span ^36 nm, 
the same distance as the average step size of myoV 
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Figure 10. Inverse pitch (|im _1 ) distributions for myoV- 
6IQand myoV-4IQ. Histograms of the inverse pitch (jim -1 ) 
calculated from individual actin filaments twirling along 
either myoV-6IQ (top) or myoV-4IQ (bottom) construct. 
The histograms are normalized. Actin filaments with an 
inverse pitch <±0.4 um -1 (pitch > 2.5 um) are classified as 
non-twirlers. myoV-6IQ twirls with a left-handed pitch, with 
a mean value of 1.4 pm, whereas myoV-4IQ twirls either 
with a left- or right-handed pitch, with mean values of 1.2 
and 1.0 um, respectively. 
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(Mehta et al., 1999). Because both motor domains of 
myoV bind stereo-specifically to subunits of actin, the 
azimuthal difference between the two heads of myoV is 
equal to the difference between the two bound actin sub- 
units. Electron microscopy (EM) studies have shown that 
a single step of myoV can span 11, 13, or 15 actin subunits 
(Walker et al., 2000; Oke et al., 2010), corresponding 
to discrete changes in the myoV azimuth of Aa L = —28°, 
0°, and +28°, respectively, values reasonably close to the 
three Aa L subpopulations for native myoV-6IQ (Aa L 
averaging -29°, -1.8°, and +33°). Thus, the polTIRF mea- 
surements indicate that the native myoV-6IQ principally 
spans 13 actin subunits and less frequently spans 11 or 15 
subunits, the same as observed in EM images obtained 
from stepping myoV (Walker et al., 2000; Oke et al., 2010) . 

The principal previously reported effect of shorten- 
ing the lever arm is that it proportionally shortens the 
step size; in the case of myoV-4IQ, the step size averaged 
24 nm (Purcell et al., 2002; Sakamoto et al., 2005). 
Based on its mean step size alone, myoV-4IQ_ should 
span, on average, approximately eight to nine actin sub- 
units, from which we would expect a Aa L distribution 
with peaks centered at Aa L = +111° and —55°. This ex- 
pectation is in contrast to our measured Aa L distribu- 
tion, which contained three subpopulations with peaks 
at Aa L = -35°, -1.4°, and +47°. What is the source of 
this discrepancy between expectation and observation? 
The answer to this most likely lies in the broad step-size 
distribution of myoV-4IQ obtained from optical trap 
measurements, which Purcell et al. (2002) noted is twice 
as wide as that of myoV-6IQ, implying that myoV-4IQ_ 
spans a larger range of actin subunits than myoV-6IQ. 
Considering that the step-size distribution measured in 
the optical trap is well populated for steps ranging from 
10 to 40 nm, this corresponds to a potential range of 
^2-13 actin subunit spans for myoV-4IQ. Consistent 
with this range, the components at Aa L = —35° and 
~1.4° we report for myoV-4IQ correspond closely to 
subunit spans of 11 and 13, also consistent with the ob- 
servations of Oke et al. (2010). In addition, we detected 
a significant number of steps of myoV-4IQ with a peak 
of Aa L = +47°, which could correspond to 2 (+28°), 4 
(+56°), or 15 subunit spans, although a span of 15 sub- 
units is less likely considering the given range of its step 
sizes (Purcell et al., 2002). Oke et al. (2010) did not 
observe two and four subunit spans by EM, which may be 
related to insufficient sampling of particles with dis- 
torted structure, as mentioned in their paper. 

It is interesting that, despite the broad step-size distri- 
bution of myoV-4IO_, it visits actin subunits within a 
narrow azimuthal range, Aa L ~ ±40°. This connection 
between myoV step-size and predicted azimuthal 
change was addressed in simulations performed by 
Vilfan (2005a), who examined likely actin-binding sites 
of myoV with lever arms containing two, four, six, or 
eight IO_ motifs. The model used in this analysis relied 



on calculating the bending energies of the four myoV 
constructs. Because of high bending energies in the 
myoV lever arm associated with large azimuthal changes, 
all of the constructs were limited to spanning 2, 11, 13, 
or 15 actin subunits, for which the azimuthal changes 
fall in the range of Aa L = ±28°. When Vilfan later in- 
cluded the natural fluctuations in the actin helix, which 
can vary ±5° in the azimuthal direction (Egelman, 
1997) , the model then predicted steps up to Aa L = ±56° 
for some of the constructs (Vilfan, 2005b) . This model 
well predicts the observed step sizes for myoV-4IQ and 
myoV-6IO_, and predicts subunit spans of 2, 4, and 1 1 for 
myoV-4IO_, similar to those implied by our Aa L distribu- 
tion (Fig. 8 B) and the myoV-4IO_ step-size distribution 
(Purcell etal., 2002). 

The mean step size (Purcell et al., 2002; Sakamoto 
et al., 2005) and actin subunit span (Oke et al., 2010) of 
myoV are clearly dependent on the length of its lever 
arm. It is therefore surprising that the means of the Aa L 
distributions, related to the span over actin subunits, 
are so similar between myoV-6IO_ and myoV-4IQ (Fig. 8 
and Table 3) . Our results on myoV-6IO_ and myoV-4IO_, 
and the Vilfan model, strongly suggest that the azi- 
muthal range of myoV on the actin filament is largely 
independent of the length of its lever arm. This con- 
straint on its azimuth is likely a result of the relatively 
high stiffness of the myoV lever arm, probably leading 
to approximately equal azimuthal stiffness as the actin, 
which would energetically exclude large bending angles 
(Vilfan, 2005a; Sun and Goldman, 2011) and may even 
straighten the actin helix somewhat between the two 
bound heads. Although the lever arm constrains the azi- 
muthal motions of myoV-4IQ to a range similar to that 
of myoV-6IQ, the myoV-4IQ construct takes more side- 
ways steps than myoV-6IQ, presumably because its trun- 
cated lever arm makes straight steps that necessitate 
stretching 36 nm, less energetically favorable. We specu- 
late that the consequent extra strain increases the likeli- 
hood of stepping sideways. Nevertheless, the range of 
azimuthal angles is limited by the discrete angles of the 
actin helix, relatively independent of the neck length. 
A different isoform, myosin VI, has a much larger range 
of step sizes and azimuthal positions because its lever 
arm is very flexible (Sun et al., 2007). 

Subpopulations in 0 P versus (p P distributions may indicate 
different labeling positions for the BR-CaM 

The distribution of 0p versus <ppfrom the processive runs 
of the myoV-6IQ constructs reveal three fairly distinct 
subpopulations (Fig. 6 A and Table 4) that quantify the 
orientation of the probe relative to the lever arm of the 
individual myoVs that were observed (Fig. 1 A). A crys- 
tal structure-based model of the isolated LCD with as- 
sociated CaMs shows that each CaM has a different 
binding orientation relative to the heavy chain (Terrak 
et al., 2005). The three different subpopulations of 
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(9 P , <p P ) pairs most likely arise from labeled CaMs bind- 
ing to different IQ motifs in the lever arm. Strikingly, 
the three subpopulations of {9 P , <p P ) in myoV-6IQ are 
also apparent in the myoV-4IQ distributions with similar 
peak values (Fig. 6 A and Table 4) . Because myoV-4IQ 
has only the first four IQ motifs in its lever arm (Materi- 
als and methods) , this similarity in the observed (9 P , q> P ) 
pairs suggests that the IQ motifs to which labeled CaMs 
are bound are limited to sites 1-4 for both constructs. 
The existence of multiple exchange sites is supported 
by other single-molecule experiments, also using BR- 
CaM-labeled myoV, from which they also concluded 
that their results came from a heterogeneous popula- 
tion of myoVs labeled at different IQ sites (Syed et al., 
2006; Topraketal., 2006). 

We can more closely identify the IQ sites that likely 
correspond to our data by using a recendy constructed 
molecular mechanical model of a myoV dimer with the 
two motor domains bound to actin 13 subunits apart 
(Parker et al., 2009). From this model, we calculated 
the parameters 9 P , <p P , (3 P , <x P , (3 L , and a L for each of the 
six CaMs in both its leading and trailing positions 
(Table S5 and Fig. SI 1) . The probe orientation was 
taken from a vector drawn between the (3-carbon atoms 
of the two amino acids in CaM replaced with cysteines 
and labeled with BR (P66 and A73) (Corrie et al., 1998; 
Forkey et al., 2003) ; the lever-arm orientation was taken 
from a vector tangent to the myosin heavy chain near 
the center of each IQ motif. The probe orientations 
corresponding to the first three IQ positions have val- 
ues for (6 P , <p P ) that are similar to the centers of the sub- 
populations from the (6 P , <p P ) distributions for myoV-6IQ 
and myoV-4IQ (Tables 4 and S5, and Fig. SI 1 ) , whereas 
CaM positions 4-6 differ substantially from the mea- 
sured (9 P , <p P ) values. This suggests that the three sub- 
populations isolated from the 9 P versus q> P distributions 
resulted from the labeled CaMs binding to the first, sec- 
ond, and third IQ motifs. The CaM bound to the sec- 
ond IQ motif of myoV has been reported to dissociate 
most easily in the presence of calcium (Trybus et al., 
2007) and so is the most likely one to exchange with 
BR-CaM in our labeling procedure (Materials and meth- 
ods). Furthermore, because the second CaM makes 
contacts with neighboring CaMs, it is possible that its 
dissociation would destabilize the CaMs bound to the 
first and third IQ motifs, thus making them targets for 
exchange with BR-CaM as well. 

Evaluating the assumptions in the myoV analysis 

Two important assumptions about myoV stepping are 
required for determining (6 P , (p P ) in the present analysis. 
The first assumption is that its lever arms are straight in 
the double-headed attached state during stepping. 
Given a myoV with two straight lever arms of the same 
length, the myosin and actin form an isosceles triangle 
and ^Ljmding + P L.traiiing = 180°. Previous studies have 



proposed that the lever arm on the leading head can 
adopt a "telemark" or kinked configuration (Snyder et al., 
2004) . However, markedly kinked lever arms are unusual 
in cryo-EM images (Walker et al., 2000; Oke et al., 2010) . 

To determine how a curved or kinked lever arm could 
perturb our analysis, we again turned to the molecular 
model of myoV (see above and Parker et al., 2009), 
which has a curved leading lever arm. The sum (5 fading + 
Putraiiing from this model can be compared with the 
180° value imposed by the analysis. In this model, the 
probe at IQ position 3 has /3 LMding + f3 UmU ,n g — 180°, and 
the values for (dp, <p P ) closely match the values deter- 
mined from the third peak of the (9 P , <p P ) distribution 
(Table S5 and Fig. Sll). p LMdlng + /3 L , (ra ,y„ !{r in IQ positions 
1 and 2 of the model are 163° and 144°, and the values 
for (9 P , <p P ) in peaks 1 and 2 of the experimental 9 P 
versus cp P distributions (Table S5 and Fig. SI 1) are also 
moderately shifted by 10° and 20°, respectively, from 
the corresponding model orientations. This difference 
illustrates that the primary consequence of the assump- 
tion of straight lever arms is a possible small deviation 
of 9 P and <p P related to the degree of curvature of the 
lever arm. 

The second assumption made in the analysis is that 
there is no azimuthal twist (7) of the lever arm around its 
own axis relative to the motor domain. 7 and <p P sum to- 
gether in all coordinate calculations, so as long as 7 is con- 
stant, its value is equivalent to setting an offset for <p P 
without loss of generality. If 7 is variable, then calculated 
cp P would be offset from its true value. Some of the offsets 
observed between the myoV model and our data (Fig. SI 1 
and Table S5) may be accounted for by twist in the lever 
arm between the leading and trailing positions. 

Extending the analysis for myoV to other systems 
The analysis described here can be applied to other bio- 
logical molecules, particularly molecular motors, rotary 
energy convertors, and nucleic acid-processing enzymes. 
Many enzymes and macromolecules undergo rotational 
motions as part of their functional mechanisms. For the 
detailed analysis described here to work, some of the 
features of molecular motor stepping should apply, in- 
cluding discrete orientations that are repeated in suc- 
cessive enzymatic cycles, a fixed or detectable rotational 
frame of reference, and labeling methodology that fixes 
a probe dipole relative to a domain that undergoes func- 
tionally relevant rotations. Other myosins, kinesins, or 
dyneins that walk along filaments with a hand-over-hand 
mechanism could use the analysis presented in this pa- 
per with little modification. The assumptions used in 
the analysis (straight lever arm and constant local twist 
angle) would need to be valid. In other myosins and in 
kinesins, the structure of the waiting state between 
steps and the dynamic changes during steps are not de- 
fined well (Block, 2007; Sweeney and Houdusse, 2010). 
With dynein, the basic motion of the step or stroke is not 
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understood well, such as whether its AAA + ring rotates 
during a step (Numata et al., 2008; Gennerich and Vale, 
2009). These applications of detailed polTIRF analysis 
are relatively straightforward. When the local orienta- 
tion of the probe (Bp, (pp) is already known (i.e., from a 
crystal structure), the assumption that the labeled do- 
main is straight could be relaxed. This was not done for 
myoV because its 12 potential labeling positions have 
different <pp angles. 

AAA + ring proteins, besides dynein, constitute a large 
family of macromolecular machines requiring function- 
ally important rotational motions. In the mitochondrial 
Fl ATP synthase, the coupling between the hydrolysis or 
reformation of ATP and tilting of the subunits is described 
by several different models (Nakanishi-Matsui et al., 2010; 
Okuno et al., 2011) that might be distinguished using 
polTIRF. Similarly, the coupling between ATPase activity 
of viral portal motors and pumping of DNA into the viral 
capsid is not setded (Moffitt et al., 2009) . 

The ribosome and ribosomal elongation factors un- 
dergo rotational motions, but the relationship between 
these motions to proofreading of aminoacyl-tRNA selec- 
tion and to translocation of the messenger RNA are 
unknown (Schmeing and Ramakrishnan, 2009) . DNA- 
and RNA-processing enzymes are natural applications 
of orientational analysis resulting from the helical na- 
ture of duplex nucleic acids. Topoisomerases and heli- 
cases exhibit large rotational motions, whose coupling 
to topological adjustments of the DNA and to unwind- 
ing of secondary structure is uncertain (Bustamante 
et al., 2011; Klostermeier, 2011). In all of these cases, de- 
tails of the structural dynamics in real time may be de- 
termined and correlated with the functional outputs by 
analyzing polTIRF recordings as we have done for myoV. 

Conclusions 

We found that the removal of two IQ motifs from the 
LCD of myoV led to myoV-4IQ taking fewer straight 
steps than the native myoV-6IQ construct, as we expected. 
However, considering the 24-nm step size of myoV-4IQ, 
the observation that the myoV-4IQ and myoV-6IQ took 
steps to the left and right with similar magnitude but 
different frequencies, as measured by Aa i; was un- 
expected. These findings suggest that the step-wise 
changes in the azimuth of myoV are limited by the stiff- 
ness of its lever arm to A.a L = ±28° and their frequency is 
determined by the length of the lever arm. We also 
found three distinct subpopulations in the distribution 
of the (0p, (pp) angles (Figs. 6 A and S9). We tentatively 
attribute these subpopulations to distinct binding sites 
of the labeled BR-CaM on the lever arm of myoV. This 
conclusion supports other studies that have also in- 
ferred BR-CaM binding to three separate sites on the 
LCD (Syed etal., 2006; Toprak etal., 2006). The (9 P , <p P ) 
values calculated from our polTIRF data also roughly 
agree with values predicted for the IQ, positions 1, 2, 



and 3 from a molecular model of myoV (Parker et al., 
2009). The analysis that we developed to transform the 
measured probe angles to the more functionally rele- 
vant lever angles is sufficiently general to be applied to 
other proteins, with repeating conformational changes 
measured using polTIRF. 

APPENDIX 

1 . Derivation of expression for the orientation of the probe 
relative to actin (a P , j3 P ) as a function of the lever-arm 
orientation (a L , /3 L ) and probe orientation relative to the 
lever-arm frame {9 Pl q> P ) 

The intensities measured in our experiments are used to 
calculate the orientation of the probe attached to myoV 
relative to the laboratory frame of reference (8 "laboratory* 
^laboratory) > which we then transform to the actin reference 
frame {a P , f3 P ). The orientation of the probe can also be 
described relative to the frame of the lever-arm (6 P , <p P ), 
whose orientation is itself expressed relative to the actin 
frame (a L , (3 L ) (Fig. 1 A and Materials and methods) . 

In deriving an expression relating the orientations of 
the probe to the lever arm in the actin frame, we note 
that this expression must also account for the orienta- 
tion of the probe relative to the frame of the lever arm. 
We construct the rotation matrix ^R^, which is derived 
from four separate rotation operations that rotate the 
probe frame (P) to the actin frame (A) (Fig. S2) . The 
combined rotation matrix A R P is derived by multiplying 
together the four rotation matrices. We describe rota- 
tions as changes in the angles (3 L , a L , Op, and (pp. The first 
two rotations that align the z A axis with the lever arm are 
performed so that the resulting axes are consistent with 
our definition of the lever-arm frame of reference (Fig. 1 A 
and Materials and methods) .These initial rotations are 
performed on the actin frame of reference, defined 
by the axes x ft y 0 , and z 0 (where x 0 = x A , y 0 = y A , and z 0 = z A ; 
Fig. S2 A) . The first rotation matrix describes the rotation 
of these axes about the z 0 axis by an amount equal to 
(it/2 — a L ) , so that the y 0 axis is aligned with the projec- 
tion of the lever arm that lies in the Xg-y 0 plane perpen- 
dicular to the actin filament: 



R =Rot\z 0 ,-\--a L 



sin a L 

-cos a, 



cos a L 0 
sin a, 0 



(A.1) 



This rotation results in a coordinate frame defined by 
x h y h and Zj, where Zj = Zn (Fig. S2 B). The second rota- 
tion matrix describes the rotation about the xj axis by 
an amount equal to fit, so that the Zj axis is aligned with 
the lever vector: 



'R 2 = Rot{x v -P L ) = 



1 0 0 

0 cos P L sin P L 
0 - sin P L cos P L 



(A.2) 
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This rotation results in a coordinate frame defined by 
x 2 , y 2 , and z 2 , where x 2 = xj = x L , y 2 = %, and z 2 = %l, where 
the Xl, yi, and Zi axes describe the lever-arm frame of 
reference (Figs. S2, B and C). The third rotation ma- 
trix describes the rotation about the z 2 axis by an 
amount equal to (tt/2 — cp P ), so that the y 2 axis is 
aligned with the projection of the probe that lies in the 
x s~y3 plane perpendicular to the lever arm: 



2 R 3 =Rot\z ,,-\ — <p f 



sin (p p cos cp p 0 
- cos cp p sin cp p 0 
0 0 1 



(A.3) 



This rotation results in a coordinate frame defined by 
X3, j3, and Zs, where Zj = z 2 (Fig. S2 D) . The final rotation 
matrix describes the rotation about the x 3 axis by an 
amount equal to 8 P , so that the Zs axis is aligned with the 
probe vector: 



3 R 4 = Rot(x 3 ,-6 p ) 

The total rotation operation is described by 



1 0 0 

0 cos 9 P sin 6 P 
0 -sin0 o cos£L 



(A.4) 



A„f O-r.4 



R 



V ■ R 2 ■ 2 R 3 



R 



(A.5) 



The result of this last rotational operation is that the 
Zs axis is aligned with the vector of the probe \ P 
(Fig. S2 E), where the X4 = x P , y 4 = y P , and Z4 = z P axes 
represent a probe frame of reference. In the probe 
frame, v P is a unit vector aligned with the Z4 axis: 



(A.6) 



v P is a vector defined in the probe frame (P) and can 
be redefined in the actin frame (A) as the vector v A : 











V P = 


v P.y 




0 











= v. 



R v, 



(A.7) 



This is converted to a^and using the following trigo- 
nometric relationships: 



a p =atan2(^,w AiJ ) 
P P = acos (^A, z ) 



(A.8) 



(A.9) 



Solving the above relationships, we find the two expres- 
sions for apand /3^as functions of a L , f3 L , Op, and (p F : 
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Pp(dp,<P P ,P L ) = acos (-sin ji L ■ sin 8 p ■ sin tp p + cos P L ■ cos 0 P ) , 

(A.10) 



where: 



0<p „<n (A.11) 



a P (6p,(p P ,a L ,P L ) = Man2(v Ax , v Aj ) , (A. 1 2) 

where expressions for v A , x and VA, y are given in the main 
text Eq. 1, and: 



-K < a „ < 71 



(A. 13) 



and the atan2(x, y) function is similar to the standard 
arctan())/x) function, except that atan2(x, y) is single 
valued over a larger range (— tt < a P < it vs. —it/2 < a P < 
it/2 for arctan (31/ x) ) . 

2. Derivation of expression for lever-arm orientation 
relative to actin (a L , /3 L ) as a function of the probe 
orientation relative to actin (a P , fS P ) and relative to the 
lever-arm (8 Pl ip P ) 

We wish to find an expression for the orientation of the 
lever arm in terms of the orientation of the probe, both 
in relation to the actin frame of reference (x A , y A , z A ). 
We start by rearranging Eq. A. 10: 



-sin fi L ■ smOp ■ sin cp p = cos P L ■ cos0 p - cos P P (A- 14) 
( - sin P L ■ sin 6 p ■ sin <p p ) 2 = ( cos P L ■ cos 9 p - cos P p ) 2 



(A.15) 



- cos P L 2 ) ■ sin8 p 2 ■ sinq> p 



= cos p, ■ cos y„ + 



cos P p - 2 ■ (cos P L ■ cos 6 P ■ cos P p ) 



(A. 16) 



2 ( 2 2 2 \ 

cos P L • I cos 0 p + sin 9 p ■ sin (p p I 

-cos P L ■ (2 ■ cos0 p ■ cos P p ) + (cos p p - s\nO p ■ sin^ 2 ) = 0 (A. 17) 
We then find an expression for /3 L : 



8, ■ cos fi p ± ^sindj 1 - sirup/ -(sin/i/ -sm9 p * -cos^') 



cos 0 + sin 8 ■ sin q> 

P P T p 



where: 



0 < P L < n. 



J 

(A.18) 



(A. 19) 



To solve for a L , we rearrange Eq. A. 12: 

tana p 

(— cos<p p + tana L ■ cos B ■ sinip^ ) ■ sinS^ + tana L ■ sin P L ■ co$9 p 
(tana L • cos(p p + cos P L ■ sin<p p ) ■ sm9 p + sin P L ■ cos9 p 

(A.20) 

tan a ■ ( sin 0 ■ tan a ■ cos <p - sin 0 ■ cos B ■ sin m - sin B ■ cos 0 ) 

L v P P 'P P " L "P "L P * 

= tana p ■ \sin$ p ■ cos fi L ■ sin^ p + sin )8 ■ cos0 p ) — sin0 p ■ tana p ■ cos(p p 

(A.21) 

Rearranging this yields an expression for 

o-l ( e p • <Pp . «p - &> ) = a tan2 ( , w A>J ) , (A.22) 

where expressions for Wa, x and Wa, } are given in the main 
text Eq. 3, where: 



is applied to a modified Eq. A.25: 

tanf * Aa p > J = tan(a f >+<1 - a p < ), 
and Eq. A. 12, simplifying to: 

tan(a f M - a P ri ) 

2 2 

tana i7i+m - sec a Lri — tana in sec a in+ 
1 — tan 2 a L n ■ tan 2 a L 

Using the trigonometric identity 



sec A = 1 + tan A, 



Eq. A.31 simplifies to: 



1 («/..„+,„ -«,.,„ ) : 



tan a , — tan a , 



1 + tan a , ■ tan a , 



(A.30) 



(A.31) 



(A.32) 



(A.33) 



-k < a, < n. 



(A.23) Using the identity given in Eq. A.29, we find that: 



3. Equality between m \a P and m \a L when m is even 

We consider the difference in orientation of the probe 

between two angular states: 



"AP n =P M -P n 



tan(a f -a p J = tan(a i ntm - a L 
which is equivalent to: 

a P,n+m ~ a p,n ~ a L,n+m ~ a i.,n 



(A.24) when the tan function is single valued; i.e.: 



(A.34) 



(A.35) 



where: 



Aa = a 



n = \,2 ...N — m 



m = 1,2 ... iV — n. 



(A.25) 



(A.26) 



(A.27) 



Note that m denotes the number of the intervals in 
the difference ™A and is even when two leading or two 
trailing states of myoV are considered. Given the condi- 
tions of our model for myoV, where 9p and cpp are con- 
stant over all n and m and the lever arms are straight 
(Materials and methods) , when m is even: 



(A.28) 



Under the conditions given above, and when m is even, 
the trigonometric identity 



i(A~B) = 



tan A - tan B 
1 + tan A ■ tan B 



(A.29) 



— <(a -<*)<—. 
2 * * * 2 



(A.36) 



We complete the derivation by noting that using Eq. 
A.25, Eq. A.35 can be identically expressed as: 



1 4 K > 

Aa p = Aa L 



(A.37) 



when m is even. 



4. Reconstructing the relative a L path using the 
relationship 2 Aa P = 2 Aa L 

We consider Eq. A.37, under the conditions given in 
Section 3 of this Appendix, when m = 2: 



When m = 1, Eq. A.25 can be expressed as: 



(A.38) 



(A.39) 



and when m = 2, Eq. A.25 is: 

2A « n =«„ + 2 -«„ =(««+! -«n) + («» + 2 -«„ + l) = A «„ + A « n+ 1- 

(A. 40) 
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In general, the difference between two angular states 
separated by m intervals can be expressed as: 



k + A J _. h h A , / = > A ■ 

n n+l n+\m—l) / , n + i 



(A.41) 



over all values of m and n. Using Eq. A.41, Eq. A.38 is 
then expressed as: 



Aa Pn + Aa Pn+1 = Aa Ln + Aa Ln+1 , 



(A.42) 



which is rearranged to define a factor, R, which ex- 
presses the effect of the probe orientation out of the 
plane of lever-arm rotation on Aa^and is constant over 
all n: 



Aa Pn - Aa L n = Aa L n+l - Aa Pn+1 



(A.43) 



= • (AOj,^, - Aa L nHm _ 1} ) = R. 

This allows us to express A.a P as the sum of A.a L and R: 
Aa Pn =(-1)" -R + Aa Ln . < A ' 44 ) 

Therefore, calculating A.a L from A.a P depends on 
finding a reliable value for R We take the average of our 
measured values of Aa P over all n: 



Aa„ 



-•(Z(-D" 



Aa p J. 



(A. 45) 



Using Eq. A. 44, this is then expressed as: 

Aa f = X((-l) Aa Pv +R) 

N -1 " 

1 A '~' 

= ^ + y<(-l)" -Aa p J. 

AT-1 ^ 
From this relation we find that: 

Aa,, =7? 

is true for all values of n and m when: 



(A.46) 



Aa, = Aa, . , 

L,n L,n + m ' 



(A.47) 



(A.48) 



and the maximum value of n = N — 1 is even, where jVis 
the total number of angular states. 

We then extend this to the general case in which Eq. 
A.48 does not hold. Because of this, the mean of these 
Ac*f6 does not necessarily equate to R as in the simple 
case detailed above. If we assume for a given run of 
myoV that Eq. A.48 is true for some of the Aa^s be- 
cause there is a predominant azimuthal motion upon 
each step, we can calculate R by taking a weighted aver- 
age of all of the Aa^ H values. Aweighted average excludes 



values of Aa ft , that significantly deviate from the aver- 
age Act/). The weighting function, W, is defined by: 



w 



(I 



Aa„ -Aa„ +e 



(A.49) 



where Aa^„ is defined in Eq. A44, and Aa p is the un- 
weighted average. In addition, a small value, e, is chosen to 
avoid division by 0, and ris an even positive integer. A value 
near zero is added to the residual in the weighting func- 
tion to prevent division by 0. The weighting function is 
then applied by combining Eq. A. 45 with Eq. A.49 to give: 



W Aa„ 



Aa rj 



(A.50) 



IT" 



Because the weighting function excludes values of 
Aaj> n that significantly deviate from the mean Aa p , and 
hence values of Aa iK that are significantly different 
from the mean Aa L (for instance, when a relatively 
straight-walking molecule takes a sudden sideways step) , 
Eq. A.48 is largely satisfied, and: 



weighted 



(A.51) 



Eq. A.50 therefore allows us to calculate a reasonable 
value for R, even in the more general case where Eq. 
A.48 does not hold for all n, as long as there are a signif- 
icant number of Aa L H s for which it is true (e.g., the mol- 
ecule either walks relatively straight [Aa L „ = 0] , or it has 
a typical azimuthal motion [Aa^ „ = &a L ]). 

5. Rationale for calculating the orientation of the probe 
relative to the lever-arm {6 Pl (p P ) for individual processive runs 

We again consider the probe at two angular states, n and 
n + m (Eqs. A.24-A.27) , in the case where the condition 
for /S^given in Eq. A.28 is satisfied. Solutions for a L and (3 L 
are defined over all 6 F -cpp space as: 



a L,n = a L{ d P><Pp' a p,„'Pp,„) 

Pl,„=Pl(8p><Pp,Pp,„) 

We first consider the general case where: 

-«L,„ = A «L,„ 



(A.52) 



(A.53) 



(A.54) 



(e p ,(p p ,a Pn ,P Pn )-a L (e p ,(p p ,a Pn+ln ,P Pll+m ) = Aa Lll 



Pt.n + Pt.n+m Pl,n 



(A.55) 



(A.56) 



P L (e P ,<P P , Pp,„ )+P L (e P ,<p P , P P ,„ +m ) = T °'P L , n , (A - 57 > 
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where m is nonzero and odd. 

We then look at the specific case where m= 1, la 'l3 L = it, 
and Aa L = 0: 

a L ( ft , <P P . «p.„ > ft,„ ) ~ a L (ft . 9? , oc P , n+1 , ft„ +1 ) = 0 (A.58) 

ft (ft ><Pp ■ ft,» ) + ft (ft ><?>p ■ ft,» + i ) = « (A M) 

Eqs. A.58 and A.59 represent a system of two equations 
that share the six parameters 0 P , cp P , (5 Kn , f3p M \, a Pn , and 
ap M i. Given our measured values for f3 Pn> (3 Pn+ i, a Pn , and 
a Pn+ i, the values for the parameters 6 P and <p P can be 
solved exacdy, assuming a real solution exists, which is ad- 
dressed below. We can also substitute A.a L „ with the values 
estimated using Eq. A.44 because the conditions under 
which this equation is used must also be met when calcu- 
lating 6 P and <p P . Because we do not have any way to di- 
recdy estimate lot ^L,m for this modified general case we 
use the assumption 7 °'/3/.„ = 180° (corresponding to the 
situation where both lever arms are straight) , as discussed 
in Materials and methods. This expression is already rep- 
resented in Eq. A.59. 

For a given set of [3 P and a P angles, a range of values of 
6 P (0 < 9 P < it) and (p P (—ir <cp P <ir) were inputinto modi- 
fied versions of Eqs. A.55 and A.59: 

(a L (e^cp^a^.Pr J-a^e^y^a^.Pr^^-Aa^ =0 

(A. 60) 

(ft (ft , <P P , ft,„ ) + ft (ft , <? r , ft,„ +I )) - x = 0 (A-61 ) 

Optimum values for 6 P and <p P corresponded to min- 
ima in plots of 9 P and <p P against the values calculated 
from Eqs. A.60 and A.61. Over the range of 6 P and ip P , 
two such minima would appear, as expected (Materials 
and methods); one set of 0pand <ppis then selected to 
calculate fi L and a L using Eqs. A. 18 and A.22. 
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